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By the method of generahzed spherical harmonic polynomials, the 
Schrodinger equation for a four-body system in D-dimensional space is re- 
duced to the generalized radial equations where only six internal variables 
are involved. The problem on separating the rotational degrees of freedom 
from the internal ones for a quantum iV-body system in D dimensions is 
generally discussed. 



I. INTRODUCTION 



Recent years have witnessed a flurry of investigations into the arbitrary D-dimensional 
problems [1-5] in many branches of physical chemistry and chemical physics. The prob- 
lems associated with the D-dimensional hydrogen atom [6-8], the D-dimensional harmonic 
oscillator [9-11], and the connection between the two [12-16] have been thoroughly dis- 
cussed. During the past few years, with the application of dimensional scaling to the 
quantum theory of atomic and molecular structure, large-D helium problem has also 
been studied by many authors [17-21]. This approach requires solving the few-body 
Schrodinger equation in a D-dimensional coordinate space and has been applied to a 
large number of physically interesting problems [22-25]. Due to the complexity of the 
problem for an A^-body system in D dimensions, so far there is no complete theoretical 
solution when > 3. 
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In our recent work [26], a new method for separating the rotational degrees of free- 
dom from the internal ones in a few-body system was proposed. The power of this new 
approach is in its ability of great simplification in calculation of energy levels of a few- 
body system in terms of the generalized radial equations involved only internal variables, 
which are derived from the Schrodinger equation without any approximation. Some typ- 
ical three-body system in three-dimensional space, such as a helium atom [27-29] and a 
positronium negative ion [30] have been solved numerically with high precision. The key 
to the approach is that we have found a complete set of independent eigenfunctions of an- 
gular momentum for the system, which are homogeneous polynomials in the components 
of Jacobi coordinate vectors and satisfy the Laplace equation, and chosen a suitable set of 
internal variables. Any wave function with a given angular momentum can be expanded 
with respect to the base functions where the coefficients, called the generalized radial 
functions, depend only upon the internal variables. The generalized radial equations sat- 
isfied by the generalized radial functions are easily derived owing to the nice property of 
the base functions [26]. This method has been generalized to the arbitrary dimensional 
space for a three-body system [31]. The exact interdimensional degeneracies in the system 
can be obtained directly from the generalized radial equations [32]. 

To further this study, we expect to apply this approach to an A^-body system in 
D dimensions. As noticed in our previous paper [31], the cases with N < D are very 
different to the cases with N > D. The general formulas are hard to express uniformly 
due to arbitrariness of D and N. However, the main characters are manifested fully in a 
four-body system of D dimensions, but not fully in a three-body system. The four-body 
problems also play a fundamental role in nuclear and hypernuclear physics [33-35]. In 
this paper we will study the problem of separating the rotational degrees of freedom from 
the internal ones for a quantum four-body system in D dimensions in some detail. The 
general case (A"-body system) will be summarized. 

The plan of this paper is as follows. In Sec. II, after separating the motion of the 
center of mass by Jacobi coordinate vectors, we will define the generalized spherical har- 
monic polynomials for a four-body system in D dimensions and prove that they constitute 
a complete set of independent base functions for a given total orbital angular momentum 
in the system. Some new features in comparison with the three-body case are also dis- 
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cussed in this section. The generahzed radial equations satisfied by the generahzed radial 
functions are established in Sec. III. In Sec. IV, we will generalize this method to sepa- 
rate the rotational degrees of freedom from the internal ones for an iV-body system in D 
dimensions. Some conclusions will be given in Sec. V. 



II. THE GENERALIZED SPHERICAL HARMONIC POLYNOMIALS 

For a quantum A'"-body system in an arbitrary D-dimensional space, we denote the 
position vectors and the masses of particles by and by m^, k — 1, 2, . . . , N , 

respectively. M = Yl,k'^k is the total mass. The Schrodinger equation for the A^-body 
system with a pair potential V , depending upon the distance of each pair of particles, 
|rj - Tfel, is 

1 ^ 

^ fe=l 
where v?^. is the Laplace operator with respect to the position vector r^. For simplicity, 

the natural units h — c— 1 are employed throughout this paper. The total orbital angular 

momentum operators Lab in D dimensions are defined as [21,36] 



N 



d d 



Lab ^ -Lba^ Vfe«^ "^kb^ — r' a,6 = 1,2, . . .D, (2) 

k=l I '^''^kb Of ka J 

where r^a denotes the ath component of the position vector r^. Now, we replace the 
position vectors with the Jacobi coordinate vectors R^: 



1<J<(A^-1), M, = ^mfe, M^ = M, (3) 

fe=i 

where Rq describes the position of the center of mass, Ri describes the mass-weighted 
separation from the second particle to the first particle. R2 describes the mass-weighted 
separation from the third particle to the center of mass of the first two particles, and so 
on. It is straightforward to illustrate that the potential is a function of Rj • R^ and is 
rotationally invariant. 

In the center-of-mass frame, Rq = 0. A straightforward calculation by replacement 
of variables shows that the Laplace operator in Eq. (1) and the total orbital angular 



momentum operator Lab in Eq. (2) are directly expressed in R^: 
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For a four-body system, there are three Jacobi coordinate vectors Ri, R2 and R3, 
which will be denoted for simphcity by x, y and z, respectively: 



X — 



777.1777.2 



1/2 



{r2 - ri} , 



z = 



1/2 



7?7i + 7772 
(7771 + 1712)1713 
777.1 + 7772 + 7773 
(7711 + 7712 + + 1714)1114 
1711+1712+ ms 



miTi + 77l2r2l 

rs \ 



(5) 
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= Vx + Vy + Vl, 
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a<b=2 



T _ r (^) I r (2/) I r (^) 

-t^ab — -t^ab + ^ab + "^afe 
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a<b=2 
D 



(6) 



a<b=2 



a<6=2 



The Schrodinger equation (1) for D > = 4 reduces to 

{ + + V^} *(x, y, 7) = -2{E-V (e„ 77,, 0)} *(x, y, z), 

6 = X • X, ^2 = ?7i = X • y, ^3 = (1 = X • z, (7) 

^2 = y • y, 7^3 = C2 = y • z, Cs = z • z. 

where ^j, rij. and arc internal variables. It is worth noticing that for the cases 3 = D < 
N two Jacobi coordinate vectors x and y can determine the body-fixed frame and this set 
of internal variables is not complete because two configurations with different directions of 
z reflecting to the plane spanned by x and y are described by the same internal variables. 
As pointed in Ref. [26], the variables (3 has to be changed to (x x y) ■ z. We will further 
discuss this problem in Sec. IV. 

Since Eq. (7) is rotational invariant, the total orbital angular momentum is conserved. 
As discussed in Ref. [31] , in D-dimensional space, the wave function ^'(x, y, z) with a given 
total angular momentum has to belong to an irreducible representation of SO{D), and the 
angular momentum is also denoted by the representation. For a four-body system, there 



are only three Jacobi coordinate vectors so that the possible irreducible representation 
is described by a three-row Young pattern [//, i^, r] of SO{D), or its highest weight M = 
(Mi,M2,M3,0,...,0), where 

Ml = ^ - i/, M2 = 1/ - T, M^^T. (8) 

We only need to consider the highest weight state ^{^[''^'^^(x, y, z) because its partners can 
be calculated from it by the lowering operators. In this paper the highest weight state 
will be simply called the wave functions with the given angular momentum [/i, i/, r] for 
simplicity. 

Now we are going to find a complete set of independent eigenfunctions of total orbital 
angular momentum, where "independent" means that each one in the set cannot be ex- 
pressed as a combination of the remaining with coefficients only depending on the internal 
variables. As discussed in our previous paper [31], the spherical harmonic polynomials 
^^^(x) are homogeneous polynomials in the components of x of degree Z, spanning an 
irreducible traceless tensor space describes by the Young pattern [l) = [Z,0, 0]. When 
D > 6, the explicit forms for some polynomials with higher weights m are as follows [37]: 

J^g(x) = Ar^,K^i + ^^2y, 

H-2,i,o,...,o) W = -VINdA^i + ix2y^\x3 + ix^), (9) 

3^(i-4,2,0,...,0) W = - ^)/^Nd,i{Xi + tX2y-\xs + tX^f, 

where the last equality holds for Z > 1, and N^^i denotes the normalization factor given in 
[37] . The product of two spherical harmonic polynomials (x) and y^) (y) belongs to 
the direct product of two representation (/) and (/'), which is a reducible representation. 

It can be reduced by the Littlewood-Richardson rule and contraction of a pair of Xa and 
Ha, where the latter relates to the internal variables: 

min{l,l'} min{l,l'}—s 

(0^(0- [l + l'-s-2t,s,0]. (10) 

s=0 t=0 

Since a base function containing a factor depending on internal variables is not indepen- 
dent, only those representations [I + I' — s, s, 0] [t — in Eq. (10)] calculated by the 
Littlewood-Richardson rule are related to the independent base functions [31]. Calculat- 
ing by the Clebsch-Gordan coefficients and removing the normalization factor, we obtain 
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the independent base functions for the representations [l+U — s, s,0], called the generalized 
spherical harmonic polynomial Qf^' '^^^{'x., y). Changing the parameters /j, — I + I' — s, 

V = s and q = I, we define the generalized spherical harmonic polynomial (^^'^(x, y) for 
the representation [p,, u] [31] as 

" {q-u)\{n-q)\ (11) 

X12 = xi + ix2, Xu ^X3 + ix4, Y12 = yi + iy2, YsA = 2/3 + «2/4- 

For the product of three spherical harmonic polynomials, Eq. (10) is generalized to 
(/)»(/')»(/") 

min{l,l'} min{(l+l'-r),(r+r)} min{r.{l" -u+r)} (12) 

~ [l + l' + r -u-T,u,r]®.... 

r=0 i'=r r=0 

The ellipsis denotes those representations related to the base functions which are not 
independent. 

Filling the digits 1, 2 and 3 arbitrarily into a given Young pattern [//, i^, r] (// > > r) 
we obtain a young tableau. A Young tableau is called standard if the digit in every 
column of the tableau increases downwards and the digit in every row does not decrease 
from left to right. In fact, the digits "1", "2" and "3" denote the components of x, y, 
and z, respectively. Obviously, the representation [I -\- V -\- V — v — r, v, r] hsted in Eq. 
(12) corresponds to a standard Young tableau, where the number of digit "1" in the first 
row is the numbers of digit "2" in the first and the second rows are respectively (/' — r) 
and r, and the numbers of digit "3" in the first, second and third rows arc respectively 
(f +r — i/ — r), u — r and r. The base functions in the remaining representation spaces, 
which correspond to non-standard Young tableaux, are not independent. 

For a given pattern [//, i/, r] , each standard Young tableau is determined by three 
parameters g, p and r, where q is the number of digit "1" in the first row, p and r are 
the numbers of digit "2" in the first and the second rows, respectively, g, p and r should 
satisfy the constraints: t < r < q and r<v<q-\-p<^. The number of standard 
Young tableaux for the given Young pattern [//, i/, r] is equal to the dimension of the 
representation [//, r] of the SU(3) group: 

d[^,.,,](5C/(3)) = !(;, - r + 2)(i. - T + l)(/i - 1/ + 1). (13) 
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For a given representation i^, r] of SO(D), each standard Young tableau denoted 
by (q, P, r) corresponds to a representation space. The highest weight state in the repre- 
sentation space {q,p,r) is the generahzed spherical harmonic polynomial (^^^^(x, y, z): 



«;(x,y,z) 



y<}—i^\^P yfJ.—g—prpr—Trpu—rrpr 
^12 -'12^12 -^12 -^13 
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{q — — g — p)!(r — r)!(z/ — r)! 

nq—rrpu- 
13 -'23 



-t^q+p— 1/ yfl—q—prpr—Trpq—rrpU—qrpT 
■'12 ^12 -'12 -'1'^ 



when q > ly, 
when q < v-i 



{q + p — uy^jjL — q — p)\{r — T)\{q — r)\{v — q)V 
T < r < q, T < v < q + p < 

X12 = xi+ ix2, Y\2 = yi + m-, = Zi+ iZ2, 

= 2^3 + ix4, Y34 = y3 + iy4, Zu = ^3 + izi, 
X5Q = + ixe, Y^G = ys + ^ye, -^56 = ^5 + ize, 

T12 — X12Y34 — X34Yi2, 7l3 = X12Z34 — X34Z12, T23 = Y12Z34 — 1^4^12, 

T = Xi2l342'56 + X34y5gZi2 + X56Yi2.^34 — Xi2y56-^34 ~ -^34^12^56 ~ -^56^34^12- 



(14) 



It is evident that Q'^^^^'x, y, z) do not contain a function of the internal variables as a 
factor, nor do their partners due to the rotational symmetry. Therefore, Q^^^.i^i., y, z) are 
independent base functions for the given angular momentum described by r] . Due 
to Eq. (12), the set of Q^^^ ^) complete. The reason why the generahzed spherical 
harmonic polynomial denoted by a non-standard Young tableau is not independent can 
be seen from the following identity: 



23^12 



7^131^2 — T^')Z 



12^12, 



(15) 



This identity is similar to the Fock's cyclic symmetry condition [38]. The left-hand-side 
of Eq. (15) corresponds to a non-standard Young tableau, and two terms in the right- 
hand-side correspond to two standard Young tableaux, respectively. 

Since the problem on completeness of the set is very important, we are going to prove 
this problem by another method. On the one hand, because the base function Q^p^l^' y-i ^) 
is a homogeneous polynomial of degree + u + t m the components of x, y and z, we 
calculate the number Rd{1) of base functions in the sets for the representations r] 
with -\-T — I. Namely, we want to calculate how many homogeneous polynomials of 
degree I exist in the sets of the independent base functions . We first calculate how many 
base functions exist in the set for a given representation r]. The dimension of the 



representation i/, r] of SO(D) is i/, r]): 

o?D([Af,i/,T]) = (L' + 2/x-2)(D + /x + z/-3)(L' + /x + r-4)(D + 2i/-4) 

X (D + + r - 5)(D + 2r - 6)(/x - r + 2)(/i - i/ + - T + 1) (16) 

{D + l2-5)\{D + u -Q)\{D + T - 7)! 
^ (L> - 2)!(D - 4)!(D - 6)\{fi + 2)!(i/ + l)!r! ' 

Thus, the number of base functions in the set for the representation u, r] is 
d[fM,u,T]{SU{3))d£){fJ',i^,T). Then, the number Ru{l) of base functions in the sets for the 
representation [p,, r] with /j, + u + t = I is: 

[1/3] [a-r)/2] 

Rd{1) = Y. E d^il-u-r),u,r]{SU{3))dD{[{l-l^-T),l^,T]), (17) 

where [x] denotes the largest integer less than or equal to x. 

On the other hand, the number of linear independent homogeneous polynomials of 
degree I in the components of x, y and z is M^Q): 

( 

Md{1) = 




After removing those polynomials in the form ^j/(x, y,z), ?7j/(x, y, z) and ('^/(x, y,z) 
where /(x, y, z) is a polynomial of degree (/ — 2), the number M^(/) reduces to Kd{1): 

Koil) = Md{1)-QMd{1-2) + 15Md{1-4)-20Md{1-Q) 
+ 15Md{1 - 8) - 6Md{1 - 10) + M^(/ - 12) 
= {{3D - 7){3D - 8){3D - 9) (3D - 10)(3L> - 11) (3^ - 12) 

+ 121{D - 4) [72 + (3D - 10) (3D - 11) (27^2 - 153D + 236)] (18) 

+ 4Z2[184 + 45(D - 4) (3D - 11) (QD^ - 57D + 98)] 

+ ASOl^{D - 4)(9D2 - 63D + 126) + 80/^(27^^ - 207D + 404) 

+ 576Z^(D-4) + 64/^}i^±^fJ|ll, 

where / + 3D > 13 and K4{0) = 1. It is checked by MATHEMATICA that 

Rd(1)^Kd{1). (19) 

Thus, we have proved again that d[fj,^i,^r]{SU{3)) polynomials (5^p^(x, y,z) construct a 
complete set of independent base functions for the angular momentum [/x, z/, rj. 
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The generalized spherical harmonic polynomial Q^!^{x.,y,z) is a homogeneous poly- 
nomial of degrees q, {p + r) and {ii + v + t — q— p — r) in the components of x, y and z, 

respectively. It is a simultaneous eigenfunction of y^, Vy) Vz) Vx ■ Vy; Vx ■ Vz, Vy ' Vz; 



and the angular momentum operators L^, 



V^g^p7(x, y, z) = v^Q^;;(x, y, z) = Vz«"(x, y, z) = 0, 

Vx ■ Vy<5^p;(x, y, z) = Vx ■ Vz<5^;;(x, y, z) = Vy • Vz(3^;;(x, y, z) = o, 

L2g^-(x, y, z) = C2([Ai, z^, T])g^;;(x, y, z), 

C2([At, i/, r]) = + L> - 2) + 1/(1/ + L> - 4) + t(t + L> - 6), 



' Q^;;(x, y, z) = g(5 + L> - 2)gp(x, y, z), 

' Q^;;(x, y, z) = (p + r)(p + r + D - 2)Q^-(x, y, z 



(20) 



(x, y,z) = {ij, + u + T - q - p - r){iJ, + v + T - q - p - r + D - 2) 
xg^-(x,y,z). 

where C2([//, i^]) is the Casimir calculated by a general formula (see (1.131) in Ref. [39]). 
The parity of Q^p^''(x, y, z) is obviously equal to (—1)''+'^+'^. 

Now, we turn to discuss the case D < 6. As is well known, the irreducible traceless 
tensor space of SO{D) described by a Young pattern has the following properties. It is 
a null space if sum of the lengths of the first two columns of the Young pattern is larger 
than D. It reduces into the self dual and antiselfdual tensor spaces if the row number of 
the Young pattern is equal to D/2. Two representations are equivalent if their Young 
patterns are the same as each other except for the first column and the sum of their row 
numbers is equal to D. Those properties cause the situation for D < 6 different to that 
for D>6. 

When D — 6, there is no problem for the representation i/, r] with t — 0, but 
when T 7^ 0, the representation is reducible. We denote the generalized spherical har- 
monic polynomials for the selfdual and antiselfdual representations by Q^p^^'^'^(x, y, z) 
and Q^p^'*'''^(x, y, z), respectively, (^^^^'''^^(x, y, z) is the same as that given in Eq. (14), 
and (5^p^'*'''^(x, y, z) can be obtained from Q^^/'''^{^,y,z) by replacing X^q, Y^q and Z56 
with X^6, and Z^g: 



^56 



- ixq, = t/s - Wq, = ^5 - ize- (21) 



The formula (16) for the dimension of the representation [/j,, v, r] of SO{D) holds for D = 6 



when T — 0. When r 7^ 0, doHfJ', v-, t\) in Eq. (16) is equal to the sum of the dimensions 
of the selfdual and antiselfdual representations such that the equahty (19) still holds for 
D = Q. 

When D = 5, in the possible Young pattern [/i, z/, r], r has to be or 1. The representa- 
tion [/X, i/, 1] is equivalent to the representation [/x, i/, 0] . Their dimensions calculated from 
Eq. (16) are also the same. The generahzed spherical harmonic polynomials Q'^'^{yi, y, z) 
given in Eq. (14) hold for = 5 except for xq — — Zq — Q and r = or 1. Therefore, 
the equality (19) holds for D = b. 

For D = 3, two Jacobi coordinate vectors, say x and y, can completely determine 
the body-fixed frame so that the variables ^3 has to be changed as (x x y) • z in order to 
distinguish two configurations with different directions of z. We have discussed in detail 
the four-body system in three dimensions in our previous paper [26]. 

The case of D = 4 is quite complicated because 50(4) is not a simple group. The repre- 
sentation [/i, z/, 0] reduces to a direct sum of a selfdual representation [{S)^^ z/, 0] and an an- 
tiselfdual one [(>l)At, v-, 0]. The generalized spherical harmonic polynomials Q^^/^'^^x.^ y, z) 
for the selfdual representations is the same as (5^^^(x, y, z) with r = given in Eq. (14). 
The generalized spherical harmonic polynomials (5^p^^'^°(x, y, z) for the antiselfdual rep- 
resentation can be obtained from Q^^l^'^^ix, y, z) by replacing X34, Y34 and Z34 with X34, 
^34 and Z'^: 

^34 = 2:3 - iX4, = y3 - iy4, ^34 = ^3 " «^4- (22) 

If r = 1, then u = 1 and the representation [/i, 1, 1] is equivalent to the representation 
[/X, 0, 0]. The standard Young tableau is described by the parameters q and p (r = 1), 
where q and p are respectively the numbers of digits " 1" and " 2" in the first row of the 
Young tableau. The generalized spherical harmonic polynomials for two representations 
[/X, A, A], A = 0, or 1, are 



^ ^ {q - \)¥if^ - q - pV-' (23) 

T = Xi2 134-^34 + -^34^34-2^12 + -^34^12-234 — -^12^34-^34 — -^34^12-^34 — -^34^34-2^12- 

The surprising thing is that Eq. (19) does not satisfied for D = 4 and I > 6. For 
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example, 

Z= 6 7 8 9 10 

i?4(0 5346 10908 20550 36332 60996 (24) 
Ki{l) 5336 10836 20256 35436 58728 

The reason is that the formula (18) for Kd{1) does not hold for D = 4 and / > 6. 

For D — A we find an identity with respect to the polynomials of degree 6 checked by 

MATHEMATICA: 

ClTla + 772^^3 + CzT^^ - 2^2Ti3T23 + 2^3^12^23 - 2773^12^13 = 0. (25) 

The identity obtained from Eq. (25) by replacing X34, ¥34 and Z34 respectively with X34, 
and Z34 still holds. Those equahties obtained by applying the lowering operators and 
(or) by multiplying a factor to above two identities are also identities. Thus, the forms 
^j/(x, y, z), ?7j/(x, y,z) and Qf{'x,y,z), where /(x, y,z) is a homogeneous polynomial 
of X, y and z of degree (/ — 2), are not independent when / > 6. It is easy to count 
by MATHEMATICA that the revised K4{1) by considering the identities coincides with 

i?4(0- 

III. GENERALIZED RADIAL EQUATIONS 

In the preceding section wc proved that ,^T-](S'f/(3)) polynomials (^^^^(x, y, z) con- 
struct a complete set of independent base functions for the angular momentum [p,, u, r] . 
Thus, any function ^jJi''^'^'(x, y, z) with angular momentum [//, i^, r] in the system can 
be expanded with respect to the base functions (^^^^(x, y, z), where the coefficients are 
functions of internal variables. 

H n—q m,in{q,v} 

*te'^'^l(x, y, Z) = ^ E ^Zri^V Vj: 0)Q^;;(X, y, Z), (26) 

q=T p=max{(i/-q),0} r=T 

where the coefficients i^ql^i^j,VjXj) are called the generahzed radial functions. When 
substituting Eq. (26) into the Schrodinger equation (5), the main calculation is to apply 
the Laplace operator (4) to the function ^|(^'^'^^(x, y, z). The calculation consists of three 
parts. In the following, we remove the arguments (x, y,z) for simplicity. 

The first is to apply the Laplace operator to the generahzed radial functions ip^^{^j, rjj, Q), 



11 



which can be calculated by replacement of variables: 

V'<;; = {^idl + An2dl + ACsdl + 2D {d^, + d,, + d^,) + (Ci + 772) dl 

+ (6 + C3) dl + (772 + Ca) 9,^3 + 46 (9^, + % + 46 (9^, + ^^3) % (27) 
+ 4773 (9,, + 9c3) 9,3 + 2773%% + 26%% + 26%%} 

where 9^ denotes d/d^ and so on. The second is to apply the Laplace operator to the 
generahzed spherical harmonic polynomials Q'^'^ ■ This part is vanishing because Q'^^ 
satisfies the Laplace equation. The last is the mixed application 

2 { (%^^;;) 2x + (%<;;) y + (%<;) z} ■ vxQp 
+ 2 { (%V'rp7) X + (%V'^;;) 2y + (%V'^;;) z} • VyQp (28) 
+ 2 { (%<;;) X + (%<;;) y + (%<;;) 2z} ■ VzQp. 
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Prom the definition (14) for Q^p^ we have 



■I qpr 



qQr 



■'qpr 



qpr 



z- Vzg^;; = (// + i^ + T-g-p-r)g^;;, 

' (p+l)(g-r) 



(/i-g-p+l)(r-T + i; 



(g-l)p(r+l)' 



when q > ly, 
{j^ - q + l)Q(gli)(p+i)r^ when g < i/, 



(g-i/ + l)Q 



flUT 



u-q 



(g+l)(p-l)r' 

V(g+l)(p-l)r 



when q > ly, 

{/jl — q — p + l){r — T + 1) 



when q < ly, 



{n-q-p+l){q-u + r -t) 

^(q-l)pr 



v-q 

(p+l)(i/-r + l) 



V(g+l)(p-2)(r+l)' 



g — z/ 



when q > ly, 



qpr 



-{'^-q + W{q-i){p+i){r-i)^ when g < i/, 
(g - i/ + l)(5(7+i)p^, when g > i/, 

(g+p-i/ + l)(g- r + l) _ {r - r + l){pL + u - 2q - p) 

^{q+l)pr 



v-q 



v-q 



X gf,+i)(p_i)(,+i), when q<v, 



y ■ S/xQqpr 



— < 



{H-q-p+ l)(5g(7-i)r + (J^ - r + l)Q^p7r-i)' ^ - ^' 

{ll-q-p+ 1)QJ7-I)r + (9 - ^ + l)<3Qp{r-l)' ^ < ^' 

(P + + (r - r + l)g^;(;+i), when g > ^, 

(g + p - i/ + l)(3j7+i)r- + - + l)<3gp(r+i)' ^hen q<v. 



(29) 



Hence, we obtain the generahzed radial equation, satisfied by the functions ipqpriC^ 0" 

V'<;; + 4g%^,7; + 4(p + r)^,,^^^; + 4(/i + z/ + r- p- g- r)d^,il^^;; 
2p(g-r + l)^ 2(/x-g-p)(r-T) ^ ,,M,.r 



+ 



+ 1 %'^(9+l){p-l)r 



g — z/ 

+ 2(g-z.)%^f;:,)(^^,)^ + 

2p(z/ — r 



g — + 1 
2(/i-y-y;)(^7-;/ + /--- + l) 



f %C-ri)(p-i)(.+i) + 2(? - + 2(/x - g - P)9,3<(;;i), 



q - u + 

+ 2(z/ - r)a^3^/'^p''(;^^) + 2p9,,3^j;"_i)^ + 2(r- - r)a^3^, 
^_2(^-l^)V;^-, for g>i/. 



flUT 

qp(r-l) 
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v'<;; + 4g%<;; + 4{p + r)d^,%;; + 4(/i + r - ^ - r)^^^^- 

+ 2(/. - 5 - p)(r - r + 1)9^3 - 2p(? - r)ae3Vr;;i)(,-i)(.+i) (30&) 
+ 2p(g - r)%V^;';_\)^, - 2(r - r)(/. - g - p + 

+ 2(r - r)a,3<;(;_,) = -2 - V) , for g = 



v'<;; + + 4(p + r)a,,V'^;; + 4(Ai + i/ + T- p- g- r)ac3V'^;; 

2(/x-g-p)(r-T) 



, 2(g + p-i/)(g-r) ^ 2(r - r)(/x + - 2g - p + 1) ^ 

+ i/-g + l ^«3V^(g-i)pr i/-g+l '^&n9-i)(p+i)('-i) 

+ 2(/. - g - v)dnA7i^.)r + 2(? - r)9,3V';7,+,) + 2(g + p - ^)9,3V'J7_,), 

+ 2(r - r)9,3^^;(;_,) = -2 (£; - y) V'^-, for g < i., 

(30c) 

where S/'^ipqpr is given in Eq. (27). Only six internal variables ^i, ^2, ^3, V2, '^3, and C3 are 
involved both in the equations and in the functions. Eq. (30) holds either for D > 6 or for 
4 < D < 6. For the latter cases some selfdual representation, antiselfdual representation, 
or equivalent representations may occur. Especially, for a four-body system in D = 4 
dimensions, the representation [//, 1,1] is equivalent to the representation [/i, 0, 0], but the 
generahzed radial equations for them are decoupled. They will be coupled for the A^-body 
system with N > D = 4. 



IV. QUANTUM iV-BODY SYSTEM IN D DIMENSIONS 

It is hard to write a unified formulas of the generalized radial equations for an A^-body 
system in arbitrary £)-dimensions. However, from the study of the three-body [26,31] 
and four-body system, we are able to summarize the main features on separating the 
rotational degrees of freedom from the internal ones for an A^-body Schrodinger equation 
in D dimensions. 

First, after removing the motion of the center of mass, there are {N — 1) Jacobi 
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coordinate vectors Rj for an A'"-body system. On the other hand, in an D-dimensional 
space it needs {D — 1) vectors to determine the body-fixed frame. When D > N, all 
Jacobi coordinate vectors arc used to determine the body-fixed frame, and all internal 
variables can be chosen as Rj ■ R^. The numbers of the rotational variables and the 
internal variables are (A^ - 1){'2D - N)/2 and N{N - l)/2, respectively. When D < N, 
only (D — l) Jacobi coordinate vectors are involved to determine the body- fixed frame, and 
the rest can be expressed by the first (D — l) Jacobi coordinate vectors and the internal 
variables. The set of internal variables Hj ■ R^ is no longer complete because it could not 
distinguish two configurations, say with different R/j refiecting to the superplane spanned 
by the first (D — l) Jacobi coordinate vectors. The correct choice for the internal variables 
are 

Cjk — Rj ■ Rfc) Ca ^ ^ ^ai...aoRlai ■ ■ ■ R{D—l)ai}-\Raai}-i 

ai...ao (31) 

l<j<^-l, J<k<N-l, D <a<N -1. 
The numbers of the rotational variables and the internal variables are D{D — l)/2 and 
D{2N - D - l)/2, respectively. 

Second, for an A'^-body system in D-dimensions {D > N), the angular momentum is 
described by an irreducible representation of SO{D) denoted by an {N — l)-row Young 
pattern \p] = [/xi, fi2, ■ ■ ■ , fJ^N-i], /^i > At2 > • • • > fJ^N-i- Due to the rotational symmetry, 
one only needs to discuss the eigenfunctions of angular momentum with the highest weight. 
The complete set of independent base functions with the highest weight consists of the 
eigenfunctions Q[g|(Ri, . . . Rat-i) identified by the standard Young tableau (q). Filling the 
digits 1, 2, . . ., N —1 arbitrarily into a given Young pattern [fi] we obtain a young tableau. 
A Young tableau is called standard if the digit in every column of the tableau increases 
downwards and the digit in every row does not decrease from left to right. Any standard 
Young tableau is described by a set of parameters (g) which contains [N — 1){N — 2)/2 
parameters g^fe, 1 < k < j < N—1, denoting the number of the digit j in the kth row in the 
standard Young tableau. The number of independent base functions Q'^|(Ri, . . . R^v-i) in 
the complete set is equal to the dimension rfj,^] [SU (A^ — 1)] of the irreducible representation 
[//] of the SU(A^ — 1) group. (5(^|(Ri, . . .Rat-i) is a homogeneous polynomial of degree 
J2 IJ-k with respect to the components oi (N — l) Jacobi coordinate vectors Rj, and satisfies 
the generahzed Laplace equations [see Eq. (20)]. The explicit form of Q^^(Ri, . . .Rd-i) 
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for the given standard Young tableau (q) is very easy to write. In the Young tableau, 
for each column with the length t, filled by digits ji < j2 < ■ ■ ■ < jt, Q[^|(Ri, . . . Rd-i) 
contains a determinant as a factor. The rth row and sth column in the determinant is 
Rjr{2s-i)+iRjr{2s) if -D > 2(A^— 1). (5[g|(Ri, . . . Rd-i) also contains a numerical coefficient 
for convenience. When N < D < 2{N — 1), some selfdual representation, antiselfdual 
representation and equivalent representations have to be considered just like the discussion 
given in the end of Sec. II. When D < N, only the first {D — 1) Jacobi coordinate vectors 
are involved in the base functions (5|^|(Ri, . . . Rd_i), which are the same as those for 
smaller N = D. 

At last, when D > N, any wave function ^|(^(Ri, . . . , Rat-i) with the given angular 
momentum [ji] can be expanded with respect to the complete and independent base 
functions (5(g) (Ri, . . . , Rat-i) 

*M(Ri,...,RAr-i)=E <!(OQB(Ri,---,Riv-i), (32) 

(«) 

where the coefficients ■0|o|(O> called the generalized radial functions, only depends upon 



the internal variables. When D < N, '0(J)(O and Q^(Ri, . . . , Rat-i) in Eq. (32) have to 



be replaced with ip\'f^{C,,C) and (5(g) (Ri, . . . ,Ri:)-i), respectively. Substituting Eq. (32) 
into the iV-body Schrodinger equation in the center-of-mass frame 

N-l 

E Vr, ^S(Ri, • • • , Rtv-i) = -2{E-V (C)} ^M(Ri, • • • , Riv-i), (33) 

one is able to obtain the generalized radial equations. The main calculation is to apply the 
Laplace operator to the function ^[(i^'^'^^ (x, y, z). The calculation consists of three parts. 
The first is to apply the Laplace operator to the generalized radial functions ijj]^^ (^) , which 
can be calculated by replacement of variables. When D > N we have 

' N-l 



N-l N-l 

+ E E [(e,. + e..) C + + ^«-) (34) 

j=l k=j+l 

N-l N-l N-l 1 

+ 2E E E <!(0, 
j=i j^k=i j^t=k+i ) 

where — Ckj and 9^ denotes d/d$, and so on. The second is to apply the Laplace op- 
erator to the generalized spherical harmonic polynomials. This part is vanishing because 
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the polynomials satisfy the Laplace equation. The last is the mixed application. When 
D > N we have 



where the formulas for • VRjQ(g) ^^'^ R-fc • VR^Q^^j can be calculated from the property 
of the polynomial (5((^|(Ri, . . . Rat-i). When D < N, the internal variables have to be 
chosen as those given in Eq. (31) so that Eq. (34) becomes more complicated and Eq. 
(35) contains more terms of ■ VRjQ(g) [26]. 



In this paper, the problem of separating the rotational degrees of freedom from the 
internal ones for the Schrodinger equation of a four-body system in D dimensions is studied 
in detail by the method of the generalized spherical harmonic polynomials. We have 
found a complete set of independent base functions with the given angular momentum 
described by an irreducible representation [//, i^, r] of SO{D). This set of base functions 
have different form for the case D > 6 and 3 < D < 6. We have provided an appropriate 
choice of internal variables for this system and derived the generalized radial equations 
depending solely on internal variables. The main features on the problem of separating 
the rotational degrees of freedom from the internal ones for the Schrodinger equation of 
a A'"-body system in D dimensions is summarized. 
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